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The quantum spin Hall (QSH) phase is a time reversal invariant electronic state with a bulk 
electronic band gap that supports the transport of charge and spin in gapless edge states. We show 
that this phase is associated with a novel Z2 topological invariant, which distinguishes it from an 
ordinary insulator. The Z2 classification, which is defined for time reversal invariant Hamiltonians, 
is analogous to the Chern number classification of the quantum Hall effect. We establish the Z2 
order of the QSH phase in the two band model of graphene and propose a generalization of the 
formalism applicable to multi band and interacting systems. 
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The classification of electronic states according to 
topological invariants is a powerful tool for understand- 
ing many body phases which have bulk energy gaps. This 
approach was pioneered by Thouless et al.|l| (TKNN), 
who identified the topological invariant for the non inter- 
acting integer quantum Hall effect. The TKNN integer, 
n, which gives the quantized Hall conductivity for each 
band a xy — ne 2 /h, is given by an integral of the Bloch 
wavefunctions over the magnetic Brillouin zone, and cor- 
responds to the first Chern class of a (7(1) principal fiber 
bundle on a torus 0,0- An equivalent formulation, gen- 
eralizable to interacting systems, is to consider the sen- 
sitivity of the many body ground state to phase twisted 
periodic boundary conditions 0,0- This topological clas- 
sification distinguishes a simple insulator from a quan- 
tum Hall state, and explains the insensitivity of the Hall 
conductivity to weak disorder and interactions. Nonzero 
TKNN integers are also intimately related to the pres- 
ence of gapless edge states on the sample boundaries 0. 

Since the Hall conductivity violates time reversal (TR) 
symmetry, the TKNN integer must vanish in a TR invari- 
ant system. Nonetheless, we have recently shown that the 
spin orbit interaction in a single plane of graphene leads 
to a TR invariant quantum spin Hall (QSH) state which 
has a bulk energy gap, and apair of gapless spin filtered 
edge states on the boundary [3. In the simplest version of 
our model (a it electron tight binding model with mirror 
symmetry about the plane) the perpendicular component 
of the spin, S z , is conserved. Our model then reduces to 
independent copies for each spin of a model introduced by 
Haldane 8], which exhibits an integer quantum Hall effect 
even though the average magnetic field is zero. When S z 
is conserved the distinction between graphene and a sim- 
ple insulator is thus easily understood. Each spin has 
an independent TKNN integer ra;. TR symmetry re- 
quires rif + ri| = 0, but the difference — ny is nonzero 
and defines a quantized spin Hall conductivity. 

This characterization breaks down when S z non con- 
serving terms are present. Such terms will inevitably 
arise due to coupling to other bands, mirror symmetry 
breaking terms, interactions or disorder. Though these 
perturbations destroy the quantization of the spin Hall 



conductance, we argued that they do not destroy the 
topological order of the QSH state because Kramers' the- 
orem prevents TR invariant perturbations from opening 
a gap at the edge[2j. Thus, even though the single defined 
TKNN number (the total Hall conductance) is zero, the 
QSH groundstate is distinguishable from that of a simple 
insulator. This suggests that there must be an additional 
topological classification for TR invariant systems. 

In this paper we clarify the topological order of the 
QSH phase and introduce a Z2 topological index that 
characterizes TR invariant systems. This classification is 
similar to the TKNN classification, and gives a simple 
test which can be applied to Bloch energy bands to dis- 
tinguish the insulator from the QSH phase. It may also 
be formulated as a sensitivity to phase twisted bound- 
ary conditions. We will begin by describing our model of 
graphene and demonstrate that the QSH phase is robust 
even when S z is not conserved. We will then analyze the 
constraints of TR invariance and derive the Z2 index. 

Consider the tight binding Hamiltonian of graphene in- 
troduced in Ref. 0, which generalizes Haldane's modelQ 
to include spin with TR invariant spin orbit interactions. 

H = tJ2 4 c i + iX so ^3 c \ sZc j (!) 
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The first term is a nearest neighbor hopping term on 
the honeycomb lattice, where we have suppressed the 
spin index on the electron operators. The second term 
is the mirror symmetric spin orbit interaction which in- 
volves spin dependent second neighbor hopping. Here 
Vij = (2/\/3)(cli x d2) 2 = ±1, where di and &2 are 
unit vectors along the two bonds the electron traverses 
going from site j to i. s z is a Pauli matrix describing 
the electron's spin. The third term is a nearest neigh- 
bor Rashba term, which explicitly violates the z — ► — z 
mirror symmetry, and will arise due to a perpendicular 
electric field or interaction with a substrate. The fourth 
term is a staggered sublattice potential (f; = ±1), which 
we include to describe the transition between the QSH 
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di t(l + 2 cos a; cosy) di2 — 2fcosssiny 

c(2 A„ di5 Aso(2sin2x — 4 sin xcosy) 

d,i \r(1 — cos x cos y) d23 — A/fCosxsiny 

c(4 — v3 Ajj sin x sin y d24 VS Xr sin x cos y 

TABLE I: The nonzero coefncents in Eq. 2 with x — k x a/2 
and y — ^/Sk y a/2. 

phase and the simple insulator. This term violates the 
symmetry under twofold rotations in the plane. 

H is diagonalized by writing </> s (R + ad) = 
M as (k)e* k R . Here s is spin and R is a bravais lattice vec- 
tor built from primitive vectors a^ 2 = (a/2)(\/3y ± x). 
a = 0, 1 is the sublattice index with d = aj//v3. 
For each k the Bloch wavefunction is a four compo- 
nent eigenvector |u(k)) of the Bloch Hamiltonian ma- 
trix Tt(k). The 16 components of 7i(k) may be writ- 
ten in terms of the identity matrix, 5 Dirac matrices T a 
and their 10 commutators T ab = [T a ,T b ]/{2i) 9\. We 
choose the following representation of the Dirac matrices: 

r (l,2,3,4,5) = f^x g j j a z ,gj 7) a y g> gX^y ^ g y ^ a y ^ g z^ 

where the Pauli matrices a k and s k represent the sub- 
lattice and spin indices. This choice organizes the ma- 
trices according to TR. The TR operator is given by 
Q\u) = i(I eg) s v )\u)*. The five Dirac matrices are even 
under TR, 0r a -1 = T a while the 10 commutators are 
odd, er^e^ 1 = -r ab . The Hamiltonian is thus 

5 5 

H(k) =^d Q (k)r a + Y, d ab(k)T ab (2) 

a=l a<b=l 

where the d(k)'s are given in Table I. Note that Ti(k + 
G) = H(k) for reciprocal lattice vectors G, so Tt(k) is 
defined on a torus. The TR invariance of Ti. is reflected in 
the symmetry (antisymmetry) of d a (d a b) under k — > — k. 

For Ar ■ — the there is an energy gap with magnitude 
|6\/3Aso-2A |. For A„ > 3%/3A S o the gap is dominated 
by A„, and the system is an insulator. 3^/3Xso > K> 
describes the QSH phase. Though the Rashba term vio- 
lates S z conservation, for A_r < 2^/?>\so there is a finite 
region of the phase diagram in Fig. 1 that is adiabatically 
connected to the QSH phase at A/j = 0. Fig. 1 shows the 
energy bands obtained by solving the lattice model in a 
zigzag strip geometry Q for representative points in the 
insulating and QSH phases. Both phases have a bulk en- 
ergy gap and edge states, but in the QSH phase the edge 
states traverse the energy gap in pairs. At the transition 
between the two phases, the energy gap closes, allowing 
the edge states to "switch partners" . 

The behavior of the edge states signals a clear differ- 
ence between the two phases. In the QSH phase for each 
energy in the bulk gap there is a single time reversed pair 
of eigenstates on each edge. Since TR symmetry prevents 
the mixing of Kramers' doublets these edge states are ro- 
bust against small perturbations. The gapless states thus 
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FIG. 1: Energy bands for a one dimensional "zigzag" strip 
in the (a) QSH phase A„ = .It and (b) the insulating phase 
A„ = At. In both cases Xso = -06t and A_r = .05t. The edge 
states on a given edge cross at ka — n. The inset shows the 
phase diagram as a function of A„ and Xr for < Xso <§C t. 

persist even if the spatial symmetry is further reduced 
(for instance by removing the C3 rotational symmetry 
in (1)). Moreover, weak disorder will not lead to local- 
ization of the edge states because single particle elastic 
backscattering is forbidden Q. 

In the insulating state the edge states do not traverse 
the gap. It is possible that for certain edge potentials the 
edge states in Fig. lb could dip below the band edge, 
reducing - or even eliminating - the edge gap. However, 
this is still distinct from the QSH phase because there will 
necessarily be an even number of Kramers pairs at each 
energy. This allows elastic backscattering, so that these 
edge states will in general be localized by weak disorder. 
The QSH phase is thus distinguished from the simple 
insulator by the number of edge state pairs modulo 2. 
Recently two dimensional versions ^3 °f the spin Hall 
insulator models [llj have been introduced, which under 
conditions of high spatial symmetry exhibit gapless edge 
states. These models, however, have an even number 
of edge state pairs. We shall see below that they are 
topologically equivalent to simple insulators. 

The QSH phase is not generally characterized by a 
quantized spin Hall conductivity. Consider the rate 
of spin accumulation at the opposite edges of a cylin- 
der of circumference L, which can be computed using 
Laughlin's argument|l2j. A weak circumferential electric 
field E can be induced by adiabatically threading mag- 
netic flux through the cylinder. When the flux increases 
by h/e each momentum eigenstate shifts by one unit: 
k — > k + 2ir/L. In the insulating state (Fig. lb) this 
has no effect, since the valence band is completely full. 
However, in the QSH state a particle-hole excitation is 
produced at the Fermi energy Ep. Since the particle and 
hole states do not have the same spin, spin accumulates 
at the edge. The rate of spin accumulation defines a spin 
Hall conductance d(S z )/dt = G s xy E, where 

G s xy = -r ((S z )l - (S z )r) \e f - (3) 
Here the expectation value of S z is evaluated for the left 
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and right moving states at Ep. Since the edge states are 
not necessarily S z eigenstates this spin Hall conductance 
is not quantized. G s xy is zero in the insulating phase, 
though, provided Ep is in the gap at the edge. If in the 
insulator the edge states cross Ep, then in a clean system 
there could be spin accumulation at the edge (resulting 
from the acceleration of the edge electrons in response to 
E). However, if the edge states are localized then there 
will be no spin accumulation. Thus the nonzero spin 
accumulation persists only for the QSH phase, justifying 
the term quantum (but not quantized) spin Hall effect. 

In the quantum Hall effect, the states with zero and one 
flux quantum threading the cylinder are distinguished by 
the charge polarization. The two states can not be con- 
nected by any operator that locally conserves charge. In 
the QSH effect there is no simple conserved quantity dis- 
tinguishing the two states. However, the states are dis- 
tinguishable, because the state with an edge particle-hole 
excitation at Ep can not be connected to the ground state 
with a local TR symmetric operator. Note, however, that 
if a second flux is added, then there will be TR invariant 
interactions which do connect the state with the zero flux 
state. This suggests that the state with one flux added 
is distinguished by a Z2 "TR polarization" . 

The classification of quantum Hall states on the 
cylinder according to Laughlin's argument is inti- 
mately related to the TKNN classification of the Bloch 
wavefunctions0. To establish the corresponding topo- 
logical classification for TR invariant systems we con- 
sider TR constraints on the set of occupied Bloch wave- 
functions |tii = i : 2 (k)) . |iti(k)) form a rank 2 vector bun- 
dle over Brillouin zone torus. TR introduces an involu- 
tion on the torus which identifies pairs of points k and 
— k. Wavefunctions at the identified points are related by 
\ui(— k)) = Q\ui(k)), implying that the bundle is "real". 
Since 6 2 = —1, 6 has period 4, so that the real bundle is 
"twisted" . These bundles are classified within the math- 
ematical framework of twisted Real K theory 0|. It is 
found that such bundles have a Z x Z2 classification on 
a torus 0|- The first integer gives the rank of the bundle 
(i.e. the number of occupied bands). The Z2 index is 
related to the mod 2 index of the real Dirac operator p^. 
In the following we will explicitly construct this Z2 index 
from the Bloch wavefunctions and show that it distin- 
guishes the QSH phase from the simple insulator. 

TR symmetry identifies two important subspaces of 
the space of Bloch Hamiltonians Tt(k) and the corre- 
sponding occupied band wavefunctions \ui(k)). The 
"even" subspace, for which QHQ^ 1 = H have wavefunc- 
tions with the property that Q\ui) is equivalent to \ui) 
up to a U (2) rotation. From Eq. 2 it is clear that in this 
subspace d a fc(k) = 0. TR symmetry requires that H(k) 
belong to the even subspace at the T point k = as well 
as the three M points shown in Fig. 2a, b. The odd sub- 
space has wavefunctions with the property that the space 
spanned by 0|u.;(k)) is orthogonal to the space spanned 



by \ui(k)). We shall see that the Z2 classification may be 
found by studying the set of k which belong to the odd 
subspace. 

The special subspaces can be identified by considering 
the matrix of overlaps, (ui(k)\Q\uj(k)) . From the prop- 
erties of 9 it is clear that this matrix is antisymmetric, 
and may be expressed in terms of a single complex num- 
ber as e,jP(k). -P(k) is in fact equal to the Pfaffian 

P(k)=Pf[(« i (k)|0|« i (k))], (4) 

which for a 2 by 2 antisymmetric matrix Aij simply picks 
out Ai2- We shall see below that the Pfaffian is the nat- 
ural generalization when there are more than two occu- 
pied bands. -P(k) is not gauge invariant. Under a U(2) 
transformation \u'j) — Uij\uj), P' = PdetU. Thus P is 
unchanged by a SU(2) rotation, but under a U(l) trans- 
formation U = e* e , P' — Pe 2lB . In the even subspace 
Q\ui) is equivalent to \ui) up to a U(2) rotation, and we 
have |P(k)| = 1. In the odd subspace P(k) = 0. 

If no spatial symmetries constrain its form, the zeros 
of P(k) are found by tuning two parameters, and gener- 
ically occur at points in the Brillouin zone. First order 
zeros occur at time reversed pairs of points ±k* with op- 
posite "vorticity" , where the phase of P(k) advances in 
opposite directions around ±k*. For A„ ^ the QSH 
phase is distinguished from the simple insulator by the 
presence of a single pair of first order zeros of P(k). The 
C3 rotational symmetry of our model constrains k* to be 
at the corner of the Brillouin zone as shown in Fig 2a. 
If the C3 symmetry is relaxed, k* can occur anywhere 
except the four symmetric points where |P(k)| = 1. The 
number of pairs of zeros is a Z2 topological invariant. 
This can be seen by noting that two pairs ±kj 2 can 
come together to annihilate each other when k* = — k^. 
However a single pair of zeros at ±k* can not annihi- 
late because they would have to meet at either T or M, 
where |P(k)| = 1. If TR symmetry is broken then the 
zeros are no longer prevented from annihilating, and the 
topological distinction of the QSH phase is lost. 

The Z2 index can thus be determined by counting the 
number of pairs of complex zeros of P. This can be ac- 
complished by evaluating the winding of the phase of 
P(k) around a loop enclosing half the Brillouin zone (de- 
fined so that k and — k are never both included). 

/ = -L I dk • V k log(P(k) + id) , (5) 
2m J c 

where C is the path shown in Fig. 2a,b. 

When A^ = (as it is in graphene) H has a C2 ro- 
tational symmetry, which when combined with TR con- 
strains the form of 7i(k), and allows P(k) to be chosen 
to be real. The zeros of P(k) then occur along lines, 
rather than at points. We find that the zeros are absent 
in the insulating phase, but enclose the M point in the 
QSH phase as shown in Fig. 2b. In this case we find 
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FIG. 2: The zeros of -P(k) in the QSH phase occur at points 
±k* for (a) A„ 7^ and on the oval for (b) A„ = 0. (c) 
|P(0, a2)| in the QSH (solid) and insulating (dashed) phases 
for a 2 x 2 supercell using parameters in Fig. 1. (d) Point 
(At, 7^ 0) and line(A„ = 0) zeros of P(a) for the 2x2 supercell. 
In (a,b,d) the solid dots are TR symmetric points, which can't 
be zeros of P, and C is the contour of integration for Eq. (5). 



that Eq. 4 also determines the Z2 index (given by 1/2 
the number of sign changes along the path C), provided 
we include the convergence factor S. Note that though 
the sign of I depends on the sign of 5, I mod 2 does not. 
We thus conclude that the QSH phase and the insulator 
are distinguished by the Z 2 index /. 

The spin Hall insulator models studied in Refs. HfHTll 
are simple insulators with 1 = 0. Their Hamiltonian, 
when expressed in the form of Eq. 2, has d a f,(k) = 0, so 
that \ui(k)) is in the even subspace and |P(k)| = 1 for 
all k. Ref. 16 introduces a model which does appear to 
exhibit a QSH effect. 

Having established the topological classification of the 
Bloch wavcfunctions we now ask whether, in analogy 
with Ref. |J the classification can be formulated in terms 
of the sensitivity of the ground state wavefunction to 
phase twisted periodic boundary conditions. Such a 
formulation will address the topological stability of the 
many body groundstate with respect to weak disorder 
and electron interactions. It also provides the appro- 
priate generalization of (4,5) for multi band Hamiltoni- 
ans. Consider a Li x L2 sample with boundary condition 
Ti + L/c, ...) = e zak '$>(..., n, ...). For concreteness we 
consider a rectangular geometry, with Li = N±(a.i + aa) 
and L2 = JVa(ai — a 2)- For non interacting electrons, 
we may view the entire sample as a large unit cell with 
N a = 4N1N2 atoms imbedded in an even larger crystal. 
Then a plays the role of k, and the occupied single par- 
ticle eigenstates 4>i(a) play the role of ttj(k). <f>i{a) form 
a rank N a bundle on the torus defined by a 1,2- The Z2 
classification can be obtained by studying the zeros of 



Fig. 2c compares |-P(<5)| in the QSH and insulating 
phases, for a 16 site sample with N\ = N2 = 2. In the 
insulating phase there are no zeros. In the QSH phase 
the structure of the zeros in Fig. 2d is similar to Fig. 
2a, b. For X v ^ the first order zeros are at points, while 
for A„ = they are on a loop. The zeros can not be at 
the four TR symmetric points. This structure persists 
in the QSH phase for any cell size. The Z2 index I can 
be computed by performing the integral analogous to (5) 
along the contour C in Fig. 2d. 

A many body formulation requires the index to be ex- 
pressed in terms of the many particle groundstate |$(a)). 
It is interesting to note that for non interacting electrons 
($(c0|e|$(5)) = det[<&(5)|e|&(S))] = P(a) 2 . This 
suggests a many body generalization 



V(*(5)|e|$(5)). 



(7) 



We suspect that with this definition the topological struc- 
ture P(a) in Fig. 2c, d will remain in the presence of weak 
electron interactions. 

To conclude, we have introduced a Z2 topological clas- 
sifiation of TR invariant systems, analogous to the TKNN 
classification of quantum Hall states. This shows that the 
QSH phase of graphene has a topological stability that 
is insensitive to weak disorder and interactions. 
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